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Uniformity of the pseudomagnetic field in strained graphene 
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We present a study on the uniformity of the pseudomagnetic field in graphene as a function of 
the relative orientation between the graphene lattice and straining directions. For this, we strained 
a regular micron-sized graphene hexagon by deforming it symmetrically by displacing three of its 
edges. By simulations, we found that the pseudomagnetic field is strongest if the strain is applied 
perpendicular to the armchair direction of graphene. For a hexagon with a side length of 1 /im, 
the pseudomagnetic field has a maximum of 1.2 T for an applied strain of 3.5% and it is uniform 
(variance < 1%) within a circle with a diameter of ~ 520 nm. This diameter is on the order 
of the typical diameter of the laser spot in a state-of-the-art confocal Raman spectroscopy setup, 
which suggests that observing the pseudomagnetic field in measurements of shifted magneto-phonon 
resonance is feasible. 


INTRODUCTION 

Graphene is well-known for its outstanding electronic 
and mechanical properties as well as their remarkable 
coupling mm- Electronically, graphene behaves as a 
high mobility semimetal [3], and mechanically, it is the 
strongest material known to mankind [4] . That graphene 
is a semimetal constraints the currently reachable on/off 
ratios for graphene transistors [5]. Therefore, much re¬ 
search is focussed on this electromechanical coupling, as 
it may be utilized to tailor the electronic properties of 
graphene. In theory, it is even possible to get the de¬ 
sired electronic properties completely via so-called strain 
engineering [i[7]. In particular, a large mechanical de¬ 
formation could induce a bandgap jSUT^. However, even 
small deformations, or equivalently, small strain fields, 
already alter the electronic properties. These small de¬ 
formations generate an effective potential for the charge 
carriers in graphene, which is similar to the vector po¬ 
tential of a real magnetic field and is therefore called the 
pseudovector potential HMS]. In analogy to a real static 
magnetic field, a pseudovector potential can generate a 
pseudomagnetic field [nHinj- This pseudomagnetic field 
is always zero, unless the applied strain is not uniform. 

A major hurdle in strain engineered graphene is 
the (experimental) determination of the pseudomagnetic 
field strength and its uniformity. Recently, it was shown 
that the pseudomagnetic field in graphene nanobubbles 
can be as strong as 300 T [20j, which is much stronger 
than one can achieve with real magnetic fields. There¬ 
fore, strained graphene offers the unique opportunity to 
study the electronic properties of graphene at extreme 
(pseudo)magnetic field strengths. Alternatively, such 
systems could allow studies to the magneto-phonon res¬ 
onance in Raman spectroscopy [21] or enable so-called 
valley-tronics [22]-[24|. Numerical work of strain fields in 
graphene is therefore of additional value, as it not only 
estimates the strength of the pseudomagnetic field, but 
also shows the uniformity of the generated pseudomag¬ 
netic field and its dependence on the lattice orientation 
with respect to the strain direction. 


In this paper, we show calculations of the pseudo¬ 
magnetic field in a regular hexagon of graphene that 
is strained simultaneously in three different directions. 
This is particularly relevant for experiments employing 
microelectromechanical systems (MEMS), which exactly 
allow for these type of strain fields [25]. As the unifor¬ 
mity of the pseudomagnetic field determines whether it 
is detectable with a local probe, e.g. a scanneling tun¬ 
neling microscope cantilever or a laser probe, we specif¬ 
ically study this as a function of the lattice orientation. 
Eor this, we consider two hexagonal geometries with side 
length L of 100 nm and 1 jam respectively. We show 
that the pseudomagnetic field is constant within a circle 
of diameter d that strongly depends on the relative ori¬ 
entation between the graphene lattice and the strain di¬ 
rection. Eor the hexagonal geometry with L = 1 jam, we 
find a diameter d of around 520 nm. In addition, we show 
that a special geometry, which was believed to generate 
a uniform pseudomagnetic field for a specific strain field 
[26] . actually does not do so. Einally, we show that an 
uniform pseudomagnetic field can indeed be generated by 
applying a linearly varying force on two opposing edges 
of a rectangle m- 

This paper is organized as follows. Eirst, we describe 
the theoretical origin of the pseudomagnetic field. Sec¬ 
ond, we discuss the methods that we used to calculate 
the pseudomagnetic field. Then, we show our results be¬ 
fore finishing with the experimental implications and the 
conclusion. 

THE PSEUDOMAGNETIC FIELD 

The electronic structure of graphene is described by a 
tight-binding Hamiltonian which only takes the interac¬ 
tion of the TT-states into account [3]: 

H = ^“ 7 "), ( 1 ) 

n=l ^ ^ 

in which k is the momentum of the electron, {dn} are the 
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FIG. 1: (a) The lattice structure of graphene, in which 0 indi¬ 
cates the rotation angle between the armchair direction and 
the x-axis. The red and blue atoms indicate the two different 
sublattices, and the dn indicates the nearest neighbor vec¬ 
tors. (b) The geometry used in our calculations is a regular 
hexagon with side length L. The force F is applied to the 
indicated edges. 

three nearest neighbor vectors (see Fig. Ef), and {tn} are 
the corresponding nearest neighbor hopping parameters. 
This Hamiltonian is a 2x2-matrix, as each unit cell of 
graphene contains two carbon atoms. As Fig. Ef shows, 
each of these atoms forms a sublattice. The energy spec¬ 
trum corresponding to Eq. 0 contains six specific elec¬ 
tron momenta k, for which the energy is zero. These six 
points are characterized by only two momenta: K and 
K', which is a consequence of the two carbon atoms per 
unit cell. 

The nearest neighbor hopping parameters {tn} change 
under a deformation of the graphene lattice, as such a 
deformation changes the nearest neighbor distance. The 
change in hopping parameters is described by 

tn = (2) 

in which a is the nearest neighbor distance without any 
deformation, to ~ 2.8 eV is the hopping parameter with¬ 
out any deformation, and /3 ~ 2 — 3.37 is the Gruneisen 
parameter [9l [28H3Q] . Please note that for the calcula¬ 
tions below, we follow ISISD] and use a Gruneisen pa¬ 
rameter oi P = 3.37. 

In case of zero strain (|dn| = a Vn), one can make 
the usual expansion around the K-point for low energy 
electrons (k = K + k, to find a Hamiltonian that is linear 
in K, 


H = hv^cr • K, (3) 

where v-p = 3atQl2h is the Fermi velocity, and cr = 
{o'x^o'y) are the Pauli matrices. 

In the case of nonzero strain, one finds a similar expres¬ 
sion for the Hamiltonian by using 6tn = tn — to. In fact, 
one finds the same Hamiltonian as in Eq. Q, but with 
an additional term that is independent of the momentum 

K\ 


( 4 ) 


in which e is the electron charge, and A = (Ax, Ay) is 
the so-called pseudovector potential, 

1 ^ 

A^ + iAy = - y . (5) 

n=l 

Equation Q immediatly explains why A is called a 
pseudovector potential, as the additional term in H mim¬ 
ics the vector potential of a real magnetic field. There¬ 
fore, we define, in analogy to a real magnetic field, a 
pseudomagnetic field: 


dAy dAx 

However, there is one remarkable difference: B changes 
sign, if one expands the Hamiltonian around the K'-point 
instead of the K-point. Consequently, the electrons in 
the two valleys feel an opposite magnetic field [19] . This 
property is important for valley-tronics since it may allow 
to make valley filters |22]-[24|. 


ROTATION OF THE LATTICE RELATIVE TO 
THE STRAIN DIRECTION 

It is a convention to define an angle 0 as the angle 
between the x-axis and the so-called armchair direction 
of graphene (see Eig. [^). If 0 = 0 and the displacements 
|u| are much smaller than a, it is well known that the 
pseudovector potential can be written in terms of the 
deformation tensor Uij [i, j = (x^y)] [T3l [3T| . 


Ax — %y) 

( 7 ) 

A (-2u ) 

(8) 

in which Uij = \(^ + 



We now rotate the lattice with respect to the x-axis 
by an angle 0 (see Eig. [^) using the standard rotation 
matrix R in two dimensions: 


. (9) 

sm(J cosU J ^ ^ 

We note that K • dn is constant under lattice rotations, 
which implies that the usual derivation of Eqs. 0 and 
Q still holds. Therefore, we only need to find the de¬ 
formation tensor for the rotated lattice [26], before we 
rotate the pseudovector potential A: 



H = hv-pcr • (k — eA) 


( 10 ) 
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and 


f ^x\ ^ f cosO sin6> \ \ 

\ Ay J — sin 6> cos 0 J \ Ay/ J ‘ 


Substituting Eqs. (10) and 0 into Eqs. 0 and 0 , 
we find an expression for the pseudovector potential for 
arbitrary lattice rotations [26] : 


“ %?/) cos(3</>) - 2uxy sm{3(f>)) (12) 
“ '^yy) sin(3())) + 2uxy cos(3^)). (13) 


Considering the pseudomagnetic field, one sees that it 
has a rotational symmetry of 120°. A rotation of the 
graphene lattice by an angle of 120° gives the same pseu¬ 
domagnetic field and a rotation of 60° effectively changes 
the sign. This implies that the pseudomagnetic field is 
zero for an angle of 30° degrees. 


METHODS 

We define a hexagonal geometry with side length L, 
as shown in Eig. [^, of which the orientation is always 
kept constant. It is only the lattice orientation of the 
graphene that is rotated with respect to the x-axis. In 
order to calculate the pseudomagnetic field we make use 
of a 2-step approach. 

In the first step we calculate the deformation of the 
geometry in COMSOL [32], for which we consider only 
in-plane strain as well as only the linear contributions to 
elasticity theory. Eor this calculation, we assumed that 
graphene has a thickness of 3.35 A, a Young’s modulus 
of 1.02 TPa, and a Poisson’s ration of 0.165 [11|33]. The 
strain is applied via a fixed displacement on the three 
edges that are indicated with arrows in Eig. Erom the 
calculated deformation field, we export the displacement 
vectors at the position of the atoms for the second step, 
in which we calculate both the pseudovector potential 
and the pseudomagnetic field. 


RESULTS 

Let us start with a hexagonal geometry of L = 100 
nm. Eigurej^ shows the strain obtained by applying a 
fixed displacement of 10 nm to the three edges for zero 
rotation angle (6> = 0). As three sides have a fixed dis¬ 
placement, the side length L of those particular sides is 
fixed to that of the undeformed geometry. In contrast, 
the side lengths of the other sides are allowed to change. 
This exactly leads to the strain pattern that we observe: 
zero strain along three edges at which we apply the fixed 
displacement and a finite strain (15%) along the other 
three edges. Eor the same reason, the strain is highest 
at the six corners of the geometry (up to 50%). Please 


note that this 50% of strain will not be reached in an ex¬ 
perimental situation, as the corners will not be so sharp 
and the boundary condition will not be so rigid. In the 
center of the geometry, the strain reduces to 0%. This is 
due to the symmetry of the applied fixed displacements, 
as the (vector) sum of them is zero. 

The panels in Eig. show the pseudomagnetic field 
for rotation angles 0 varying from 0 to 30 degrees. A pos¬ 
itive (negative) pseudomagnetic field is indicated in red 
(blue). We observe two contributions to the total pseudo¬ 
magnetic field: the one at the six corners of the hexagon 
and the one from the center of the geometry. The pseu¬ 
domagnetic field at the corners does not change when 
rotating the graphene lattice. It remains constant with 
values around -1-60 T and — 60 T. The origin of this pseu¬ 
domagnetic field is our boundary condition (see above). 
As a consequence, we get a large relative displacement 
at the six corners of the geometry, which results in the 
observed pseudomagnetic field. 

In contrast, the pseudomagnetic field in the center of 
the geometry does depend on the rotation angle 0. The 
maximum pseudomagnetic field 5max in the center of the 
geometry is around 40 T for ^ = 0° and decreases rapidly 
when increasing the rotation angle 6>. At a rotation angle 
of 0 = 30°, the pseudomagnetic field in the center of 
the geometry is even zero, as the atomic displacement is 
parallel to the nearest neighbor vectors. Eigurej^ shows 
the pseudomagmetic field in the center of the geometry 
normalized by Bmax as a function of the rotation angle 0 
as well as the normalized pseudomagnetic field that one 
obtains from Eqs. (12) and (13). As Eig. shows, our 
numerical calculation is in perfect agreement with the 
theoretical prediction. 

The uniformity of the pseudomagnetic field is of cru¬ 
cial importance for the possible utilization of it in exper¬ 
iments and devices. To quantify the uniformity of the 
pseudomagnetic field 5, we calculate the average value 
of 5 in a circle of diameter d that is centered in the ge¬ 
ometry (see Eig. 0). The result is shown in the upper 
panel of Eig. for various rotation angles 0. In order 
to compare the uniformity of the pseudomagnetic field 
at different rotation angles 6>, we normalize each curve 
with respect to the pseudomagnetic field Bq ai d = 2 nm 
(see the lower panel of Eig. |^). Please note that the 
curves for different 0 are very similar and only deviate 
for large d. Einally, Eig. 0 shows the standard devia¬ 
tion of the normalized pseudomagnetic field cfbibq as a 
function of d. To characterize the homogeneous area of 
the pseudomagnetic field one has to set an upper limit 
on the standard deviation. We choose a standard devi¬ 
ation smaller than 1% of the maximum pseudomagnetic 
field in the center. By using this definition, we obtain a 
diameter of uniformity of around 34 nm for ^ = 0. 

Eigurej^ shows an overview of a feasible MEMS actua¬ 
tor that allows to experimentally strain graphene sheets. 
The arms of the MEMS actuator are placed under an 
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angle of 120 degrees and are connected with a rectangu¬ 
lar sheet of graphene, which has a width of 5 /iui and a 
length of 7 fim. This configuration generates a force field, 
which has the same threefold symmetry as for the hexag¬ 
onal structure above (see Fig.[^). However, the effective 
side length is now ^ 1 jam. The strain in the graphene 
sheet resulting from the same fixed displacement as be¬ 
fore is shown in Fig. ii. Although the strain in the center 
region of the geometry is on average 3.5%, the maximum 
strain value is 8% on the corners of the hexagon and 5.5% 
on the sides. The maximum strain in the entire geom¬ 
etry occurs at the positions where the graphene/MEMS 
device boundary has a corner. At these locations, the 
strain is even around 20%, which is the amount of strain 
at which the graphene sheet is supposed to start to rup¬ 
ture HIM]. Nevertheless, it is expected that such high 
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strain values will not be reached in an experiment. 

The pseudomagnetic field for different rotation angles 
of the lattice with respect to the strain directions is cal¬ 
culated and depicted in Fig. [^. The pseudomagnetic 
field has the same characteristics as the pseudomagnetic 
field of the strained hexagon in Fig. The strength of 
the pseudomagnetic field is highest in each corner of the 
hexagon and the field in the center is strongest for ^ = 0. 
Figure shows the strength of the pseudomagnetic field 
as a function of the diameter d of a circle that is placed 
identically to that in Fig. ii. For d = 0, the pseu¬ 
domagnetic field is approximately 1.2 T in the center. 
If we compare this with the maximum pseudomagnetic 
field in the hexagon with side length L = 100 nm, we 
see that it is approximately 30 times weaker. This is due 
to the following two effects. Firstly, by increasing the 
side length L by a factor of 10, the gradient in the pseu¬ 
dovector potentials is 10 times smaller, which results in 
a 10 times smaller pseudomagnetic field. Secondly, the 
applied strain is now approximately ^3.5%, which is ap¬ 
proximately 3 times smaller than the 10% that we had 
before. Combining both, explains the ^ 30 times weaker 
pseudomagnetic field. 



FIG. 2: (a) The strain in the geometry for zero rotation an¬ 
gle 0. There are three edges with zero strain and three with 
15% of strain. The corners of the geometry exhibit the largest 
amount of strain (up to 50%), whereas the strain reduces to 
0 % in the center, (b) The pseudomagnetic field as a func¬ 
tion of the rotation angle 0. At 0 = 30° the pseudomagnetic 
field disappears in the center of the geometry, as the atomic 
displacement is parallel to the nearest neighbor vectors, (c) 
The pseudomagnetic field in the center of the geometry is in 
excellent agreement with that from Eq. (12) and (p^. Please 


note that the pseudomagnetic field is normalized with respect 
to the field at ^ = 0°. 


FIG. 3: (a) The average pseudomagnetic field and its corre¬ 
sponding normalized standard deviation are calculated within 
a circle of diameter d that is centered in the geometry, (b) 
The upper panel shows the pseudomagnetic field as a function 
of the diameter d for different rotation angles. The bottom 
panel shows that all these curves fall almost on top of each 
other if one normalizes each curve with respect to the value 
of the pseudomagnetic field in the center, (c) The standard 
deviation of the pseudomagnetic field as a function of the di¬ 
ameter d for the rotation angles shown in (b). The generated 
pseudomagnetic field is uniform within d = 34 nm for 0 = 0°. 
The inset shows a zoom for small d. 
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If we now consider the uniformity of the pseudomag- 
netic field in such MEMS devices, as shown in Fig. |^, 
we obtain a diameter of uniformity of around 520 nm 
for ^ = 0°. If we would choose a standard deviation 
smaller than 5% of the maximum pseudomagnetic field 
in the center of the geometry, the diameter of uniformity 
increases to ^ 1100 nm for 6> = 0°. As both values are 
comparable to the typical laser spot size in state-of-the- 
art confocal Raman spectroscopy experiments, it may be 
possible to quantify the pseudomagnetic field strength in 
an experiment as it induces a shift in the MPR. 


Apart from the experimentally accessible geometries 
that we discussed so far, certain geometries with cor¬ 
responding force fields have been suggested that should 
give rise to a completely uniform pseudomagnetic field. 
In particular, Guinea et al m give two of such (ana¬ 
lytic) geometries with corresponding force fields. Using 
our approach, we confirmed the uniform pseudomagnetic 
field in the circular disk for the displacement field given 
in Ref. [26]. Let us next focus on the second geometry 
given in Ref. m- This geometry should give an uniform 
pseudomagnetic field when only perpendicular forces are 
applied to its boundary. The analytic expression for this 
geometry in polar coordinates is derived under the as¬ 
sumption of pure shear stresses (see the Supplementary 
Information of [26]): 



FIG. 4: (a) The computed strain field in a graphene sheet that 
can be reached with a state-of-the-art MEMS device shows 
values (~ 20%) close to the maximum strain that graphene 
can withstand. The force F is applied to the indicated edges, 
(b) The corresponding pseudomagnetic field for various rota¬ 
tion angles 0. (c) shows the pseudomagnetic field as a function 
of d and (d) shows the corresponding standard deviation. In 
this case, the pseudomagnetic field is uniform within d = 520 
nm. 


r{4>) = 


c 

[(cos 0/2 ^ sin 0/2) (±1 + 2sin0)]^^^ 


(14) 


In this equation, r is the radius, 0 is the polar angle, and 
c is a constant. The ± sign indicates two equivalent ge¬ 
ometries that are rotated by 60 degrees with respect to 
each other. It does not matter which sign we take, the 
radius always becomes imaginary for specific angles. For 
simplicity, we choose the + sign for our geometry. Fig¬ 
ure df shows both the real part (black) and imaginary 
part (red) of the radius r as a function of the angle 0 
for c = 1 /im. Please realize that r(0) is imaginary and, 
therefore, physically not allowed. We can obtain the ge¬ 
ometry shown in [26] by taking the absolute value of r(0) 
(see Fig. [^). However, if we calculate the pseudomag¬ 
netic field for this given geometry and force field using 
our numerical approach, we observe that it is not uniform 
(see Fig. [^). However, in the arms of this geometry, the 
pseudomagnetic field looks almost uniform. Based on 
this, it can be shown that a rectangle with linearly vary¬ 
ing forces perpendicular to two opposing boundaries (see 
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FIG. 5: The solution for the radius r for c = 1 /im as a func¬ 
tion of the angle 0 (a) and its corresponding geometry for \r\ 
(b). In contrast to what was thought, the resulting pseudo¬ 
magnetic field is not uniform (c). (d) The linearly varying 
forces (left panel) on two opposing boundaries of a rectan¬ 
gle that generate an uniform pseudomagnetic field. The right 
panel shows the corresponding deformation of the geometry. 
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left panel in Fig. [^) does give rise to an uniform pseu- 
domagnetic field [27]. The corresponding deformation of 
the geometry is shown in the right panel of Fig. [^. 

EXPERIMENTAL IMPLICATIONS 

A particular experiment, in which the strength of a 
tunable pseudomagnetic field may be directly probed, is 
the shift of the magneto-phonon resonance (MPR) in Ra¬ 
man spectroscopy [351438] . In order to exclude other ef¬ 
fects, the pseudomagnetic field must be uniform within 
the laser spot size in the confocal Raman spectroscopy 
experiment, which is exactly what we showed and dis¬ 
cussed in this paper. By following the recent experiment 
of Neumann et al. [21], where MPRs at low magnetic 
fields have been observed, we can estimate the exper¬ 
imental requirements for observing MPR shifts thanks 
to a tunable pseudomagnetic field. The most prominant 
MPR below 20 T is the so-called Ti transistion at around 
3.7 T, which results in an well resolved peak of the Ra¬ 
man G-line width as function of magnetic field. This res¬ 
onance peak has an experimentally observed peak line- 
width (FWHM) of around 700 mT. Assuming that we 
need for unambigously proving the presence of a pseudo¬ 
magnetic field at least a peak shift of half the line-width 
of this magneto-phonon resonance. This results in a re¬ 
quired pseudomagnetic field of around 350 mT, which 
corresponds to a strain value of 1% in our proposed de¬ 
vice. These values are by a factor 3.5 below the values 
discussed above (see Fig. [^) allowing for a non-perfect 
lattice alignment and making this approach promising for 
investigating tunable pseudomagentic fields. 

CONCLUSION 

We studied the strength and uniformity of the pseu¬ 
domagnetic field in a hexgonal sheet of graphene as a 
function of the rotation angle 0 between the armchair di¬ 
rection of graphene and the strain direction. We showed 
that the relative orientation of the graphene lattice with 
respect to the strain direction is extremely important, 
as the pseudomagnetic field even disappears for spe¬ 
cific angles. The pseudomagnetic field is strongest when 
the strain is parallel to the armchair direction in the 
graphene. To characterize the homogeneous area of the 
pseudomagnetic field, we set an upper limit on the stan¬ 
dard deviation in the pseudomagnetic field of 1%. Us¬ 
ing this definition, the pseudomagnetic field is constant 
within a diameter of ^ 520 nm for a hexagonal graphene 
sheet with a side length of 1 fim. Our results show 
that the pseudomagnetic field is detectable with a lo¬ 
cal probe such as a scanneling tunneling microscope can¬ 
tilever but even also with a typical confocal laser probe. 
One may even think of observing a shift in the magneto¬ 


phonon resonance due to the strain-induced pseudomag¬ 
netic field. In addition, tunable pseudomagnetic fields 
are important for valley-tronics since they may allow to 
make valley filters. 


ACKNOWLEDGEMENT 

We acknowledge funding from the ERC (GA-Nr. 
280140) and the EU Elagship-Graphene (contract no. 
NEGT-IGT-604391). 


[1] M. Huang, T.A. Pascal, H. Kim, W.A. Goddard, and 
J.R. Greer, Nano Lett. 11, 1241 (2011) 

[2] M. Trif, P. Upadhyaya, and Y. Tserkovnyak, Phys. Rev. 
B 88, 245423 (2013) 

[3] M.I. Katsnelson, Graphene: Carbon in Two Dimensions 
(Cambridge University Press, New York, 2012) 

[4] C. Lee, X. Wei, J.W. Kysar, and J. Hone, Science 321, 
385 (2008) 

[5] F. Xia, D.B. Farmer, Y-M. Lin, and P. Avouris, Nano 
Lett. 10, 715 (2010) 

[6] V.M. Pereira and A.H. Castro Neto, Phys. Rev. Lett. 
103, 046801 (2009) 

[7] T. Low, and F. Guinea, Nano Lett. 10, 3551 (2010) 

[8] Y. Hasegawa, R. Konno, H. Nakano, and M. Kohmoto, 
Phys. Rev. B 74, 033413 (2006) 

[9] V.M. Pereira, A.H. Castro Neto, and N.M.R. Peres, 
Phys. Rev. B 80, 045401 (2009) 

[10] F.M.D. Pellegrino, G.G.N. Angilella, and R. Pucci, Phys. 
Rev. B 81, 035411 (2010) 

[11] G. Cocco, E. Cadelano, and L. Colombo, Phys. Rev. B 
81, 241412(R) (2010) 

[12] S-M. Choi, S-H. Jhi, and Y-W. Son, Phys. Rev. B 81, 
081407(R) (2010) 

[13] H. Suzuura and T. Ando, Phys. Rev. B 65, 235412 (2002) 

[14] K-I. Sasaki, Y. Kawazoe, and R. Saito, Prog. Theor. 
Phys. 113, 463 (2005) 

[15] M.I. Katsnelson, and K.S. Novoselov, Solid State Com- 
mun. 143, 3 (2007) 

[16] F. de Juan, J.L. Manes, Marfa A.H. Vozmediano, Phys. 
Rev. B 87, 165131 (2013) 

[17] M.R. Masir, D. Moldovan, and F.M. Peelers, Solid State 
Commun. 175-176, 76 (2013) 

[18] N. Abedpour, R. Asgari, and F. Guinea, Phys. Rev. B 
84, 115437 (2011) 

[19] K-J. Kim, Ya.M. Blanter, and K-H. Ahn, Phys. Rev. B 
84, 081401(R) (2011) 

[20] N. Levy, S.A. Burke, K.L. Meaker, M. Panlasigui, A. 
Zettl, F. Guinea, A.H. Castro Neto, and M. F. Crommie, 
Science 329, 544 (2010) 

[21] C. Neumann, S. Reichardt, M. Drogeler, K. Watan- 
abe, T. Taniguchi, B. Beschoten, S.V. Rotkin, and C. 
Stampfer, Nano Lett. 15, 1547 (2015) 

[22] A. Rycerz, J. Tworzydlo, and C.W.J. Beenakker, Nature 
Phys. 3, 172 (2007) 

[23] M.H.D. Guimaraes, A. Veligura, P.J. Zomer, T. Maassen, 
I.J. Vera-Marun, N. Tombros, and B.J. van Wees, Nano 
Lett. 2012, 3512 (2012) 



7 


[24] Y. Jiang, T. Low, K. Chang, M.L Katsnelson, and F. 
Guinea, Phys. Rev. Lett. 110 , 046601 (2013) 

[25] H.H.P. Garza, E.W. Kievit, G.F. Schneider, and U. 
Staufer, Nano Lett. 14, 4107 (2014) 

[26] F. Guinea, M.L Katsnelson, and A.K. Geim, Nature 
Phys. 6, 30 (2010) 

[27] F. Guinea, A.K. Geim, M.L Katsnelson, and K.S. 
Novoselov, Phys. Rev. B 81, 035408 (2010) 

[28] A.J. Heeger, S. Kivelson, J.R. Schrieffer, and W-P. Su, 
Rev. Mod. Phys. 60, 781 (1988) 

[29] M.A.H. Vozmediano, M.L Katsnelson, and F. Guinea, 
Phys. Rep. 496, 109 (2010) 

[30] A.H. Castro-Neto and F. Guinea, Phys. Rev. B 75, 
045404 (2007) 

[31] J. L. Manes, Phys. Rev. B 16, 045430 (2007) 

[32] http://www.comsol.com 


[33] O.L. Blakslee, D.G. Proctor, E.J. Seldin, G.B. Spence, 
and T. Weng, J. of Appl. Phys. 41, 3373 (1970) 

[34] K.S. Kim, Y. Zhao, H. Jang, S.Y. Lee, J.M. Kim, K.S. 
Kim, J-H. Ahn, P. Kim, J-Y. Choi, and B.H. Hong, Na¬ 
ture 457, 706 (2009) 

[35] M.O. Goerbig, J.-N. Euchs, K. Kechedzhi, and V.I. 
Eal’ko, Phys. Rev. Lett. 99, 087402 (2007) 

[36] C. Eaugeras, M. Amado, P. Kossacki, M. Orlita, M. 
Kuhne, A.A.L. Nicolet, Yu.I. Latyshev, and M. Potem- 
ski, Phys. Rev. Lett. 107, 036807 (2011) 

[37] S. Goler, J. Yan, V. Pellegrini, and A. Pinczuk, Solid 
State Commun. 152, 1289 (2012) 

[38] C. Eaugeras, P. Kossacki, A.A.L. Nicolet, M. Orlita, M. 
Potemski, A. Mahmood, and D.M. Basko, New J. Phys. 
14 095007 (2012) 


